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Abstract 

Power-law singularities and critical exponents in n-vector models are considered 
from different theoretical points of view. It includes a theoretical approach called the 
GFD (grouping of Feynman diagrams) theory, as well as the perturbative renormal- 
ization group (RG) treatment. A non-perturbative proof concerning corrections to 
scaling in the two-point correlation function of the ip^ model is provided, showing 
that predictions of the GFD theory rather than those of the perturbative RG theory 
can be correct. Critical exponents determined from highly accurate experimental data 
very close to the A-transition point in liquid helium, as well as the Goldstone mode 
singularities in n-vector spin models, evaluated from Monte Carlo simulation results, 
are discussed with an aim to test the theoretical predictions. Our analysis shows that 
in both cases the data can be well interpreted within the GFD theory. 
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1 Introduction 

Critical phenomena in interacting many-particle systems are associated with cooperative 
fluctuations of a large number of microscopic degrees of freedom. Various physical quan- 
tities typically exhibit power-law singularities in vicinity of the phase transition point. 
These are described by the critical exponents, which are exactly and rigorously known for 
a class of exactly solved models [U 121 [3] . These are mainly the two-dimensional lattice 
models. For three-dimensional systems, exact results are difficult to obtain, and approx- 
imate methods are usually used - see, e. g., [H O El [71 |8| for a review of the applied 
here standard perturbative renormalization group (RG) methods. A general discussion 
of aggregation processes and critical phenomena in many-particle systems can be found, 
e. g., in [9l[T0]. Recently, the role of quantum fluctuations in critical phenomena has been 
considered [llj . For a general review, one has to mention that phase transitions described 
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by power laws and critical exponents can be observed in variety of systems, such as social, 
economical, biological systems, as well as vehicular traffic flow [12 \ \13 \ [T ^ I15j. 

An alternative theoretical approach in determination of the critical exponents in the 
ip'^ model has been proposed in [16]. We will call it the GFD theory, as it is based on 
certain grouping of Feynman diagrams. Moreover, a qualitative analysis is performed 
here |16j without cutting the perturbation series. More recently, this approach has been 
generalized [iTj to study the (p^ model below the critical temperature, where the so called 
Goldstone mode singularity (see, e. g., [ISl [H EQl [2ll [Ml IS] ) is observed. It refers 
to certain models, further called n~vector models, in which the order parameter is an 
n~component vector with n > 1. This provides one more example of power law behavior, 
exhibited by the transverse and longitudinal correlation functions in the ordered phase. 
Moreover, according to the recent Monte Carlo (MC) simulation results \25\ [26 l [27], it is 
very plausible that this behavior is described by nontrivial exponents, as predicted in [17j . 
We will discuss this question in Sec. [5l including some new Monte Carlo simulation results. 



2 Critical exponents of the GFD theory 

Here we discuss the critical exponents of the GFD (grouping of Feynman diagrams) the- 
ory [16], since these are important in our further tests and discussions. We consider a ip"^ 
model with the Hamiltonian 



H 



I {rop\^) + c(Vvp(x))2 + n/(x)) , (1) 



where the order parameter '/^(x) is an n-component vector with components V3i(x), de- 
pending on the coordinate x, T is the temperature, and ks is the Boltzmann constant. 
As in ([E]) and ([IZ]), it is assumed here that the order-parameter-field (pi{x) does not 
contain the Fourier components V3j(k) with A; > A, i. e., there exists the upper cut-off pa- 
rameter A. It is well known that the ip'^ model with certain dimensionality n of the order 
parameter belongs to the same universality class as the n-vector lattice spin model with 
spins being n-component vectors of unit length (see Sec. [5|). Therefore, the considered 
here predictions for the critical exponents refer also to these n-vector models. 

It has been claimed in |16j that possible values of the exact critical exponents for the 
ip^ model in two {d = 2) and three {d = 3) dimensions are 

^ d + 2j + 4m _ ^ ^ 2(1 + m)+j 
^ d(l + m + j) - 2j ' d{l + m + j) - 2j ' ^ ^ 

Here 7 is the susceptibility exponent, v is the correlation-length exponent, m may have a 
natural value starting with 1, and j is an integer equal or larger than —m. Other critical 
exponents can be calculated from these ones, using the known scaling relations. These 
values agree with the known exact solutions of the two-dimensional Ising model (m = 3, 
j = 0) and of the spherical model {j/m — )• 00). Based on the idea that d might be 
considered as a continuous parameter in ([2|) within 2 < d < 4, a prediction has been made 
also for the three-dimensional (3-D) Ising model: 7 = 5/4 and v = 2/3, corresponding 
to m = 3 and j = 0, as in the two-dimensional case. This value of is consistent with 
the logarithmic singularity of specific heat (according to a = 2 — du = 0) proposed earlier 
by Tseskis [28j. The exponents 7 = 5/4 and v = 2/3 have been later conjectured for the 
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3-D Ising model by Zhang [29]. The disagreement of these exponents with those of the 
perturbative RG method can be understood based on a critical analysis [30]. As explained 
in |16] . the equations ([2]) are meaningful for positive integer n, and we can have j = j{n) 
and m = m{n) in the case where the order parameter is an n-component vector. 

A relevant question is how ^ can be related to specific models, i. e., one needs to find 
j{n) and m{n). It is well known that the spherical model corresponds to the limit n — )• oo, 
so that we have j{n)/m{n) — >■ cx) at n — >■ cxd to obtain the known critical exponents of the 
spherical model in this limit. Consider now the critical exponent 

7 4 - d 

71 = 2-- = — ^ , 3 

' 1/ 2(1 + m) + j ' ^ ^ 

which describes the ~ k^'^^^ singularity (at A: ^ 0) of the Fourier-transformed two-point 
correlation function at the critical point. Since j{n) /■m{n) — )• oo, the asymptotic behavior 
of r/ = ry(n) for large n is 

4 - d 

i']{n) ~ — — — - at n ^ oo . (4) 

A general conjecture for j{n) and m(n) can be made assuming that the critical ex- 
ponents are given by analytic expressions valid for each positive integer n. It results in 
an essential restriction on possible asymptotic form of r/(n). In fact, the only reasonable 
possibility is r]{n) oc at n — )• cxd, where s must be a positive integer, taking into ac- 
count that j{n) has an integer value at any positive integer n. The latter means also that 
j[n) an'' holds at n — >■ oo with an integer coefficient a. An appropriate conjecture is 
i{n) = Vs{n) and m{n) = Qr-(n), where Vs{n) and Qr{n) are polynomials of orders s and 
r < s with integer coefficients, since m(n) and j{n) are integers at any integer n > 1 and 
j{n) /m{n) — )• oo holds at n — )■ oo. If we accept a physically reasonable idea that critical 
exponents and their derivatives with respect to n behave smoothly and monotonously at 
n — > oo, like some power of n, then this conjecture is the only possible one. Hence, we 
have nr](n) = T{x), where J-'{x) is some analytic function of the argument x = 1/n, which 
can be expanded as T{x) = ao + aix + a2X^ + . . . at small positive x. Here ao 7^ holds 
if s = 1, whereas s = 2 is a special case of oq = and ai 7^ 0. Similarly, s = 3 is a special 
case of ao = ai = and 02 7^ 0, and so on. Since we currently do not see any reason why 
the coefficient oq has to be zero, we consider the choice s = 1 as the most natural one, and 
we will further test only this possibility. Thus, we have j{n) = Vi{n) = j(l) + a (n — 1) 
and m(n) = QQ{n) = m(l). Here j(l) = and m(l) = 3 hold according to the prediction 
for the Ising model. It yields m = 3 and j{n) = a(n — 1) with some integer coefficient 
a. By definition, j{n) > —m holds for n > 1 and, therefore, a > 1. Only at a = 1 we 
have a bijective relation between j and n, i. e., each nonnegative integer j corresponds to 
one positive integer n and vice versa. From this point of view, the choice a = 1 seems to 
be the most natural one, since in other cases a question arises why only each second (at 
a = 2), each third (at a = 3), etc., integer value of j has a meaning. We will further test 
just the possibility a = 1, i. e., the conjecture 

m = 3 ; j{n) = n — 1 . (5) 

Not only the critical exponents of leading singularities, but also corrections to scaling 
are important in our tests. According to [16j, corrections to scaling can be represented 
by an expansion of correction factor (amplitude) in integer powers of t^i^-i ^^d t^"<~'^'^ at 
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t — )• 0, where t is the reduced temperature. Since (27 — dv') /{2iy — 7) is an integer number 
according to and 2 — = r] holds according to the known scahng relation, we obtain 
the expansion in powers of t^'^ . Thus, the correction-to-scaling exponent is 

9 = r]v = ujv (6) 

if the first expansion coefficient is nonzero. Here uj = 0/v is the correction-to-scaling 
exponent, describing the ~ corrections to the critical correlation function. It shows up 
also in the finite-size scaling analysis. Allowing that some of the expansion coefficients are 
zero, we can have = iriv, where i = i{n) is a positive integer. A conjecture i{n) = 3 + n 
together with ([5]) has been tested earlier in [31], based on the experimental data near 
A-transition point in liquid helium. This conjecture is based on the idea that 6 for some 
quantities tends to 1 at d — > 2 in the Ising case of n = 1, as well as at d = 3 in the limit 
of the spherical model n ^ 00. 

A smooth crossover to the 2Z)-Ising-behavior at n = 1 and d ^ 2, or to the spherical- 
model-behavior at n — 7- 00 can be expected in the ip'^ model. However, it is not necessary 
for the existence of such a crossover that the correction-to-scaling exponent tends to the 
corresponding limit value. A simple possibility is that the expansion coefficients for the 
nontrivial correction terms vanish in these limit cases. Moreover, even this is not the 
necessary condition at n ^ 00. In fact, if we consider a correction factor of the form 
1 + Y^fLi 6£(n)t''(")''(")^ where M is an arbitrarily large constant, then it is sufficient that 
lim„_j.oo ^^j^ I b£{n) |< C holds for any M > 1, where C is a finite constant, which 
is independent of M. It ensures that this correction factor tends to a finite constant at 
n ^ 00 for any given nonzero value of t, since r]{n) — > in this limit. It means that the 
coefficients of the included here nontrivial correction terms do not need vanish at n — )• 00. 
As regards the n = 1 case at d = 2, it is further shown in Sec. [3] that nontrivial correction 
terms can be present in the if"^ model and absent in the 2D Ising model. 

According to this discussion, Eq. ([6|) represents the most natural conjecture for the 
correction-to-scaling exponents 9 and uj, since in this case we do not need to assume that 
some of correction terms always vanish. This conjecture is supported by the recent Monte 
Carlo simulation results ^32j for the 3D Ising model on very large lattices with linear sizes 
up to L = 1536. Besides, according to the numerical transfer matrix calculations in [33], 
a nontrivial correction to finite-size scaling with the exponent uj = tj, probably, exists 
in the two-point correlation function even in the 2D Ising model. The relation ([6|) is 
incompatible with the predictions of the perturbative RG theory. However, an essential 
non-perturbative proof will be provided in Sec. [3l showing that corrections to scaling of 
the GFD theory rather than those of the perturbative RG theory can be correct. In Sec. U 
we will test how well the conjectures ([5]) and ^ are consistent with the experimental data 
very close to the A-transition point in liquid helium. It refers to the case n = 2. 

3 Corrections— to— scaling theorem and a non-perturbative 
proof 

Consider now the model 1^ at T > Tc (i. e., tq < roc), where Tc is the critical 
temperature. We consider the case where rg is the only parameter which depends on T, 
and the dependence is linear. In the following, we use the known thermodynamic relations 
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for the free energy F, specific heat Cy, internal energy U and entropy S at a fixed volume 
V: 

as well as the well known relation F = —ksTln Z, allowing to determine the free energy F 
from the microscopic model by calculating the partition function Z = J exp[—H/ {kBT)]Dip, 
where the symbol T>(p indicates that the integration takes place over all allowed configu- 
rations of (^(x). It yields 



fc<A 



where G(k) = (| V3i(k) p) (for any i = 1, 2, . . . , n) is the Fourier-transformed (according 
to (/7j(x) = V^^^"^ ^^f^^^ipiCk) e^^^) two-point correlation function. 

Let us denote by and Cy^^ the leading singular parts of U /V and Cy/V in the 

thermodynamic limit V — )• oo, represented in terms of the correlation length ^ at ^ oo, 
i. e., at T — 7- T^. From dH) and ([9]), where rg is a linear function of T, we obtain 

(/. \ sing / r \ si'^g 

/ G(k)dk) =(/ [G(k) -G*(k)](ik) , (10) 
Jk<A J \Jk<A J 

where (•)**"'3 generally denotes the leading singular part of the corresponding quantity in 
brackets, and G*(k) is the correlation function at T = Tc- The non-singular constant part 
G*(k)(ik is subtracted for convenience. 
Further on, we assume that Cy has either the usually expected power-law singularity 
or, more generally, a power-law singularity with logarithmic correction of 
the form C^'"^ oc {InO^^"^"- Since ^ ~ t-" holds at t ^ 0, Eq. dZD then leads to 
C^"''^ oc ^y'^U''''3_ As it is well known, the critical Ion g-wave fluctuations are responsible 
for the leading singularities near T^. Hence, U^^'^^ cannot be altered by a short-wave 
contribution to (llOp . so that 

(p \ sing 

/ [G(k) -G*(k)]rfk ) (11) 
Jk<A' J 

holds for any finite A' < A, i. e., the leading singularity is independent of the upper 
integration limit, which is formally set to A'. Note, however, that G(k) and G*(k) always 
correspond to the true upper cut-off A. Summarizing these relations, we obtain 

/ p \ sing 

C^ng^^l/ul / [G(k) -G*(k)](ik . (12) 

\Jk<A' J 

Since only the small-A; contribution is relevant, it might be well justified to use the 
scaling hypothesis for G(k) and G*(k), which is valid for small k and large Namely, we 
have 

G(k) = J^e^^-^^V^'sK^a, (13) 



5 



where ge{k^) are scahng functions, = holds and the term with 1 = describes the 
leading singularity, whereas the terms with i > 1 represent other singular contributions 
with correction exponents Oi > 0. Here we do not include possible analytic correction 
terms, since they cannot give the leading singularities of Cy and U. The critical correlation 
function 

G*(k) = ^6^A;(-T+^^)/'' (14) 

i>0 

is obtained at ^ ^ oo, so that we have 

(7^(z) ~ at z^oo. (15) 

For complete formal correctness, one should note that the use of (|13p and (|14p is justified 
if it yields a A'-independent Cy^^ , since in this case the calculated Cy^^ is not modified 
by the short-wave (not-small A;) contribution, which is evaluated only approximately. 
Now we are ready to formulate the main result of this section as the following theorem. 



Theorem. If the leading singular part of specific heat Cy"'^ in the actually considered 
99^ model has the form Cy""^ oc {h\^)^^°'/'^ (with A = corresponding to the usual power- 
law singularity), if this singularity is provided by the A'-independent small-k contribution 
to with the scaling hypothesis (Eqs. hlS^) and {l-j^ ) being valid for G'(k) and G*(k), 
and if ^ + 1 — a — dv > holds, then the two-point correlation function contains a 
correction-to-scaling term with certain exponent 

0^ = j + l-a-du, (16) 

corresponding to one of the terms with i > 1 in il3\}. 

Proof. Inserting ()13p and ()14p into ()12p and changing the integration variable to 7/ = /c^, 
we obtain 

(\ sing 

where 

z 

Fiiz) = I y''-%{y)dy with Uv) = 9iiy) " b,y^~^+'^y'' . (18) 


Although we are interested only in the leading singularity of specific heat, the correction 
terms have to be retained in p7p at this step, since some contributions can vanish after the 
integration and, therefore, a term with i > can be important. Let us first assume that 
the leading singular part of Cy/V is provided by a single term with certain i. In this case, 
C^^"^ has the required form C^'"^ oc (InC)^^"/'^ only if Fe{A'^) ~ [ln(A'e)]^(A'0'' holds at 
^ — 7- 00 with some exponent /i. Note that A = /i = is always true if the integral in (jlSp is 
convergent at z = 00. Such giiy), which ensures this property, certainly exists. It can be, 
e. g., any function which decays as y~'^ at ?/ ^ 00 with a > d. Obviously, the leading term 
in [ln(A',f)]'^(A'^)'' = [In A' + ln,f]'^(A',^)'^ is A'-independent only if (U = holds, in which 
case we obtain Cy^^ oc (In^)''*^^'^^^^^'''"^^)/'^. It is consistent with the required form at 
the condition (fT6|) . Such ge{y), which gives A 7^ and fi = also exists: it can be any 
function, decaying as (Iny)^^^y^'^ at y — )■ 00. The possibility ^ = is excluded, since (fT6|) 



6 



and 7 + 1 — a — dz^>0 cannot be simultaneously satisfied at £ = 0, as = holds by 
definition. Thus, the leading asymptotic term of the correlation function gives vanishing 
contribution to Cy"^. It means that either lim2_>.oo Fo(^) = holds, or the term with 
£ = gives an analytic (e. g., constant or ~ t) contribution. 

Let us now consider a possibility that two different terms with indices £ and i' are 
equally important, i. e., proportional to each other at — )■ oo, and Cy""^ oc 
holds. It is possible only if ^^(A'^) ~ [ln(A'^)]^(A'0'^ and F^'(A'^) ~ [ln(A'0]^(A'^)^' 
hold at ^ ^ oo with [i — Qnjv = fj.' — Og'/i^, implying that ^ n' (since 6i 7^ 6i' holds 
by definition). However, the obtained result for Cy""^ is A'-dependent, as it contains a 
factor of the form A"^ + BA"^ with B ^ 0. Similarly, the result is A'-dependent if any 
larger number of terms are equally important. Consequently, only one term with certain 
£ > 1 contributes to Cy^^ at the conditions of the theorem. Hence, (|13p contains the 
corresponding correction-to-scaling term with 61 given by (|16|) . □ 

The condition ^ + 1 — a — dv > is very meaningful, as it reduces to 7 > 1 according 
to the well known hyper-scaling hypothesis 

a + du = 2 , (19) 

and 7 > 1 really (or almost surely) holds for the 99^ model within 2 < d < 4. Hence, our 
theorem has some important consequences listed below. 

1. If the hyper-scaling hypothesis (fT9]l holds, then the statement (fTHl) reduces to 

0^ = 7 - 1 . (20) 

2. Since 9e is one of the correction exponents, we have 9<9i = 'y + l — a — du for the 
leading correction-to-scaling exponent 9. It reduces to < 7 — 1 according to the 
hyper-scaling hypothesis, this statement being valid if 7 > 1. 

3. The actual consideration allows a possibility that Cy""^ oc (In^)''*^"/'^ holds with an 
arbitrary value of A. However, if the decay of gii{y) at y ^ 00 is power-like, i. e., 
y~" , then we have either A = (at o" > d) or A = 1 (at a = d). 

4. Since the ip'^ model belongs to the Ising university class at n = 1, we have 9i = 
7 — 1 = 3/4 in two dimensions at n = 1, according to the known exact result 7 = 7/4 
of the 2D Ising model. Hence, a non-trivial correction to scaling with 9£ = 3/4 exists 
(if the conditions of the theorem are satisfied) in the correlation function of the (p^ 
model at n = 1 and d = 2. According to the known exact results (see, e. g. [34j ) . 
such a correction does not appear in the correlation function of the 2D Ising model 
on an infinite lattice. Apparently, the 2D Ising model is a special case, where the 
non-trivial corrections to scaling usually (but, probably, not always |33j ) vanish. 

5. The statement that the term with £ = gives vanishing contribution to Cy"'^ implies 
the existence of some cancellation mechanism in ()12p . 

The second consequence that ^ < 7— 1 hods at 7 > 1 is inconsistent with the predictions 
of the perturbative RG theory, whereas the corrections to scaling of the GFD theory, 
discussed in Sec. [21 completely agree with the proven here theorem. There are no doubts 
that the conditions of this theorem are very reasonable from the physical, as well as 
mathematical, point of view. Hence, the corrections to scaling of the GFD theory rather 
than those of the perturbative RG theory can be correct. 
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4 Best experimental evidences for the power— law singular- 
ities near the critical point 

There are a lot of different experimental evidences for power-law singularities near phase 
transition points. However, it is not our aim to give an exhaustive review of this topic, so 
that we will focus only on the best experimental evidences available. These are basically 
the specific heat measurements in zero-gravity (space) conditions [35] very close to the 
A-transition point in liquid helium. These measurements are done with a high degree 
of accuracy much closer to the critical point than in any other experiments or numerical 
simulations. Due to this reason, it is widely accepted to consider them as crucial tests of 
validity of the theoretical predictions for the critical exponents. We will briefly discuss 
also the second-sound velocity measurements of superfluid fraction [361 EZ] near the A- 
transition point in liquid helium. It is widely accepted that the A-transition is described 
by the critical exponents of the n-vector model with n = 2. 

It has been found in [35] that the experimental specific heat (Cp) data for a wide range 
of reduced temperatures 5 ■ 10"^'' < i < 10~^ below the A-transition temperature Tx can 
be well fit to appropriate ansatz of the perturbative RG theory, providing an estimate 
a = —0.0127 lb 0.0003 of the specific heat critical exponent, in a satisfactory agreement 
with the usual RG values. An alternative ansatz, including logarithmic correction and 
critical exponents a = —1/13 and 9 = 5/13 {6 = {n + 3)ryi^ within the GFD theory) has 
been proposed in [31]. Here we will test a different possibility: a = 2 — du = —1/13 and 
d = 1/13 given by ([2]), ([5]) and dSD at n = 2. 

We start our analysis with a critical reconsideration of the fits, obtained by assuming 
the usual RG correction-to-scaling exponent 9 ~ 0.529, i. e., the same one used in [35j . 
where two slightly different ansatz 



have been considered (in somewhat different notations than here) with t = 1 — T/T\ being 
the reduced temperature at T < T^. In our fits, we have used the raw data of [35], as 
well as certain binning of these data, described in (31]. This binning procedure differs 
slightly fro that one used in [35j. In fact, the binned data correspond within the error 

{l+b)t 

bars to Cp{t} = (bt)-^ J Cp{T)dT with t = t{l + 6/2) and some constant b, defining 



the averaging interval of the binning procedure. It is easy to verify that, if Cp{t) is given 
either by (f2T]) or by (p2|) . then Cp{t) also is described by the corresponding ansatz with 
only slightly different coefficients A, ai and a2- Therefore, the raw data and also the 
binned data can be fit to determine the critical exponent a. Since the variation of Cp 
within one binning interval is very small, the raw data and the binned data lie practically 
on the same curve, and the fit results are consistent within the error bars. We have used 
both data sets to verify the robustness of our fitting procedures and related analyses. 



The fit of raw data to (I2ip over the whole range of the reduced temperatures 4.75 x 
10-1° < t < 9.52 X 10-3 with fixed d = 0.529 yields a = -0.01263(20), Cq = 460.5(6.2), 
A = —447.6(6.1), ai = —0.0156(13) and 02 = 0.3306(96), where the standard errors are 




(21) 



and 




(22) 
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Table 1: The estimates ai and Q2 of the critical exponent q, obtained by fitting the Cp 
data (for ai) and the Cp = Cp — A [ait^ + a2t^^) data (for 02) to the ansatz ([23|) within 

^ ^ tmax ■ 



tmax 


Oil 


a2 


3.18 X 10"^ 
1.03 X 10^2 
3.70 X 10"^ 
1.07 X 10"^ 
2.95 X 10"^ 
1.01 X 10-5 


-0.013950(58) 
-0.012371(87) 
-0.01194(17) 
-0.01216(30) 
-0.01135(53) 
-0.00924(99) 


-0.012578(58) 
-0.012724(87) 
-0.01270(17) 
-0.01277(30) 
-0.01179(53) 
-0.00955(99) 



indicated in brackets (meaning a = -0.01263 ± 0.00020, Co = 460.5 ± 6.2, etc.). These 
values are well consistent with those reported in [35j. 

We have performed certain test of validity of such a fit, based on the following idea: 
this fit result for a should be consistent with the one obtained from a simpler ansatz 

Cp = Co + At-"^ (23) 

by fitting the data within t < tmax at so small values of tmax, at which the correc- 
tions to scaling ait^ and a2t'^^ become negligible. The effect of these correction terms 
is evaluated as ai — 02, where ai and 02 are obtained by fitting the Cp data and the 
Cp = Cp — A [ait^ + a2t'^^) data, respectively, to the ansatz (j23p . The latter data are 
obtained by subtracting correction terms, evaluated from the overall fit to ()2ip . As soon 
as the difference between ai and 02 becomes much smaller than the error bars, a good 
agreement of ai and Q2 with the value a = —0.01263(20) of the overall fit is expected, if 
this fitting procedure with 6 = 0.529 is valid. The fit results for ai and 02 depending on 
tmax are represented in Tab. [H As we can see, the values of ai and 02 very well agree 
within the error bars for the smallest tmax value 1.01 x 10"^. A good agreement with the 
estimate a = —0.01263(20), however, is not observed at this t. The disagreement cannot 
be reasonably explained by a possible inaccuracy in the Tx value. The discrepancy be- 
comes smaller only by 0.0005 if we shift T\ by —0.5 nK within the experimental error bars. 
Moreover, the shift T\ ^ Tx — 0.5 nK makes the deviations of the binned data points from 
the overall fit curve quite remarkable for 10 smallest t values (see Fig. [T]). Namely, such 
a shift increases the average deviation for these data points from —1.10(70) percents to 
—1.73(70) percents and, therefore, is not well justified. Hence, the actual fitting procedure 
with e = 0.529 is doubtful. 

In fact, if we assume that 6 ~ 0.529 really holds, then the fit to simple ansatz (j23p 
at very small t values, such as t < 1.01 x 10"^, should be considered as a more reliable 
method than the overall fit within t < 9.52 x 10^'^ with two correction terms included. 
Indeed, it is expected that these corrections, as well as higher-order correction terms are 
negligible in the first case, whereas this is not surely true for the neglected higher-order 
corrections in ()2ip at t ~ 10~^ used there. From this point of view, a = —0.00924(99) is 
almost the best estimate obtained by us, assuming 6 ~ 0.529. It does not well agree with 
the perturbative RG values, e. g., a = —0.01294 it 0.0006 reported in [38] . 
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Figure 1: The percent deviation of the measured (binned) Cp data points from the fit 
curve, obtained by fitting these data to ([2T]) with fixed exponents a = —1/13, 6 = 1/13 
(solid circles) and a = -0.01284, 6 = 0.529 (x). 



The same test can be performed for the ansatz (122p . yielding not better results: we 
obtain a = —0.01320(21) from the overall fit to (j22p . this estimate being not well consistent 
with ai = -0.00924(99) and 02 = -0.00970(99) obtained from 1^ at tmax = 1-01 x 10"^ 
These tests have been performed also for the fits of the binned Cp data, giving similar 
results. For example, the overall fit to (j2ip in this case yields a = —0.01284(34), whereas 
the fits to ([23]) give us ai = -0.0081(13) and 02 = -0.0085(13) at tmax = 1-09 x 10"^ 
The bin-averaging has been performed with equal weights for all data points of one bin, 
which is the reason for slightly larger statistical errors than for the raw-data fits. 

A reasonable explanation of the failure in the above tests is provided by the non- 
perturbative analysis in Sec. [3l Namely, the correction-to-scaling exponent 9 is likely to 
be not larger than 7 — 1. The latter quantity does not exceed 0.32 according to both the 
perturbative RG estimate 7 = 1.3169 ±0.0020 of [39j and our result in Sec.[2]7 = 17/13 ~ 
1.3077 for n = 2. Thus 9 should be remarkably smaller than 0.529, which means that the 
influence of the correction terms is still not negligible at t about 10~^. 

In the following, we test how well the data are described by the ansatz (j2ip with 
a = —1/13 and 6* = 1/13 proposed by ([2]), ([5]) and Here we consider deviations from 
fit curves, therefore the binned data are appropriate for the analysis, as they are less noisy 
than the raw data. We have found that the data can be well fit with fixed a = —1/13 
and 9 = 1/13 within t < 5.71 x 10~^, as it is evident from Fig. [H where the percent 
deviations from the fit curve are shown by solid circles. The fit is less perfect for two 
largest t values, which are omitted here. However, since 9 = 1/13 is quite small, it can 
be well explained by an infiuence of higher-order correction terms. For comparison, the 
percent deviations in the case of fixed exponents a = —0.01284 and 9 = 0.529 considered 
before are shown by symbols "x". In fact, the percent deviations are practically the same 
in both cases for t > 10"^. The average deviation for 10 smallest t values is —0.10(70) at 
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Table 2: The parameters of the binned-data fits to (12 ip within t < tmax with 9 = —a. 





a 


Co 


A 


ai 


0.2 


6.99 


-0.069(17) 


201(31) 


-474(29) 


-1.19(22) 


0.66(23) 


5.60 


-0.069(18) 


201(33) 


-474(30) 


-1.19(23) 


0.66(25) 


3.93 


-0.080(20) 


185(26) 


-473(19) 


-1.36(33) 


0.86(40) 


2.76 


-0.069(30) 


200(51) 


-461(30) 


-1.15(46) 


0.62(46) 


1.74 


-0.091(27) 


175(27) 


-494(44) 


-1.60(57) 


1.18(80) 



a = —1/13 and 6 = 1/13, i. e., remarkably smaller in magnitude than —1.10(70), obtained 
at a = —0.01284 and = 0.529. Thus, it might be true that the t — )• asymptotic is 
better described by the GFD exponents than by the perturbative RG exponents. 

We have considered a series of fits of the binned Cp data within t < tmax with fixed 
correction-to-scaling exponent 9 = 1/13 and a as a fit parameter. Since 6 = —a holds 
at n = 2 according to ([21), (l5|) and ([6]), we have performed also fits with 6 = —a. These 
fits give similar results as those with 9 = 1/13, but the statistical errors are remarkably 
smaller for very small tmax values. Therefore, they provide a more precise test of the 
GFD theory. We have found that such fits are sufficiently (acceptably) stable within 
1.74 X 10~^ < tmax < 6.99 X 10-^ as shown in Tab. El The fitted a values in Tab. [2] 
are not very accurate. Nevertheless, they perfectly agree within the error bars with the 
expected theoretical (GFD) asymptotic value —1/13 ~ —0.07692. The fits at three largest 
tmax values in Tab. [2] are sufficiently stable with respect to the temperature shift T\ — >■ 
TxzbO.S nK, and the shifted a estimates agree within the error bars with —1/13. We have 
fitted also the raw data to (j21|) with 9 = —a. Although the results are somewhat less stable 
in this case, the agreement for very small tmax values is observed, e. g., a = —0.054(27) 
at tmax = 3.93 X 10-^ a = -0.071(24) at tmax = 2.20 x 10"^ and a = -0.071(33) at 

tmax = 1.59 X 10-5. 

In fact, at so small 9 as 1/13, we practically cannot find such a tmax value at which the 
corrections to scaling are negligible. Even including two correction terms, the ansatz is 
expected to be really accurate only at very small t values. Hence, the agreement with the 
theoretical prediction a = —1/13 observed here at rather small t values is an evidence in 
favor of the result ([2]) of the GFD theory and our conjectures ([5]) and ([6]). The latter one 
might be more appropriate than 9 = (n + 3)r/i/ used in [3T|, since now it is not necessary 
to assume the existence of a logarithmic correction for a good fit of the data. 

Consider now the superfluid-fraction data of |36l |37j . It has been shown in [31j that 
the effective exponent feg, extracted from these data and plotted versus 

^5/13 

mcreases 

above the usual RG value ^ 0.67 and apparently converges to a value near 9/13 (the 
GFD theoretical value at n = 2), if fitted to a parabola. It is an evidence in favor of 
the conjecture 9 = {n + 3)7/1/, the exponents ry and u being consistent with ([2]) and ([5]). 
This behavior can be explained assuming ([6]), as well. Namely, first few coefficients in 
expansion of the superfluid fraction in powers of can be very small, in which case the 
scaling of the data in a wide range of t values is described by an effective correction- 
to-scaling exponent 9es ~ i'qi' with i > 1. Taking into account the experimental errors 
and uncertainty in the critical temperature (it depends on pressure, which is not strictly 
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constant in non-zero gravity conditions), I can be, e. g., 5, 4, or even 3 within the error 
bars. Owing to the recent MC analysis of the 3D Ising model (n = 1) on very large 
lattices [32], the conjecture ^ is more plausible than 6 = {n + ^)rjv. 

In summary, it is possible to give a reasonable interpretation of the actually discussed 
experimental and numerical (MC) data, assuming that ([2]), ([SJ and ^ hold. 



5 Goldstone mode singularities in the 0(n) models 

Here we consider a class of n-vector spin models, where the spin is an n-component 
unit vector with n > 2. These are also called 0(n) models due to the 0{n) global 
rotational symmetry exhibited by such n~vector model in absence of the external field. 
The Hamiltonian % in this case reads 




where Sj is the spin variable of the zth lattice site, /3 is the coupling constant, and h is 
the external field. The summation takes place over all pairs {ij) of the nearest neighbors 
in the lattice. Each spin has the longitudinal component s\\^ which is its projection on 
the external field, and n — 1 transverse components, which are perpendicular to the field. 
We are interested in the magnetization per spin M = (sy), as well as in longitudinal 
and transverse correlation functions. In the coordinate representation, the longitudinal 
(G|j(x)) and the transverse (G_l(x)) correlation functions are defined by 

G||(X2-X1) = (s||(xi)s||(x2))-M2 (25) 

^^(xa-xi) = (s±(xi)sx(x2)) , (26) 

where s± is any one of the transverse components. Due to the symmetry of the model, 
the correlation functions depend only on the coordinate difference X2 — xi . The Fourier- 
transformed longitudinal and transverse correlation functions are 



G||(k) = iV-i^G||(x)e-*''^ (27) 

X 

G^(k) = iV-i^G^(x)e-^''^. (28) 

X 

Consider now the behavior of an 0(n) model below the critical temperature, i. e., at 
(3 > Pc, in the thermodynamic limit L — )• oo. In this case, the magnetization M{h) and 
the correlation functions exhibit Goldstone mode power-law singularities: 

M{h) - M(+0) (xhP at /i ^ , (29) 
G_L(k) oc fc-^-L at /i = +0 and , (30) 

G||(k) cx A;~^ll at /i = +0 and ^ . (31) 

According to the standard theory [20l EH [22l [23l [M], A_l = 2 and Ay = 4 - d hold for 
2 < d < 4, and p = 1/2 is true in three dimensions. More nontrivial universal values are 
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expected according to [T7], such that 

(i/2 < A_L < 2 , (32) 

A|| = 2A_L - d , (33) 

P = {d/X±) - 1 (34) 

hold for 2 < d < 4. These relations have been obtained in |17j by analyzing self-consistent 
diagram equations for the correlation functions without cutting the perturbation series. 
Apart from the mathematical analysis, reasonable physical arguments also have been 
provided there to show that Aj^ = 2 could not be the correct result for the XY model 
within 2 < d < 4. 

The relations ()32p and (j33p are confirmed by MC simulation results for the longitudinal 
and transverse correlation functions in the 3D 0(4) model [27], where an estimate Aj_ = 
1.955 ±0.020 has been found. It has been stated [27] that the behavior of the longitudinal 
correlation function is well consistent with Ay about 0.9, in agreement with (I33p at A_l 
about 1.95, but not with the standard-theoretical prediction Ay = 1. According to (j34p . 
we have 1/2 < p < 1 in three dimensions. A reasonable numerical evidence for this relation 
has been obtained in [26] from the susceptibility data of the 3D XY model, providing the 
MC estimate p = 0.555(17) for the n = 2 case. It corresponds to Aj_ = 1.929(21) according 
to dSl). 

The longitudinal and transverse correlation functions for the 3D XY model have not 
been simulated in the mentioned here papers. Here we fill this gap, providing a new MC 
evidence that Ay < 1 holds. The simulations for several linear lattice sizes L < 512 at 
h = 0.000875, h = 0.0004375 and h = 0.00021875 have been performed at /3 = 0.55 (the 
critical coupling being /3c ~ 0.4542 [40.]), and the correlation functions (in (100) direction) 
have been evaluated, following the method described in [27] . As in [27], we evaluate the 
effective longitudinal exponents Aeff(A;) from the linear log-log fits within [A:,2/i;]. The 
effective exponents depend also on h and L, therefore plots of Aeff(A;) for different values 
of these parameters are compared in Fig. [2] to judge about the asymptotic exponent Ay, 
corresponding to the limit lim lim lim . The finite-size effects increase with decreasing 

of h. We observe that the plots of the effective exponent for L = 512 and L = 384 lie 
almost on top of each other at the smallest h value, and the plot for L = 256 also shows 
only small deviations. Thus, the thermodynamic limit is practically reached at L = 512 
for the actual values of the field h. Furthermore, the finite-/i effects are quite small in the 
considered range of wave vector magnitudes k > k2o, where ki = 2TTi/L with L = 512. 
Namely, the plots at two smallest h values almost agree within the statistical error bars. 

If we assume the corrections to scaling proposed by the standard theory, then the 
asymptotic exponent Ay can be quite well evaluated by fitting these smallest-/i data as 
functions of k. Recall that these corrections are represented by an expansion in powers 
of k'^~'^ and k'^~'^ [I8l |T9], corresponding to the expansion Aefr(A;) = Ay + X]j>i^i^"' 
three dimensions. Accordingly, we have plotted Aefr(A;) vs k in Fig. [2] and have fit these 
plots by a polynomial of k. The cubic fits well describe the shape of the theses curves 
within k € [A:20i^7o] and yield Ay = 0.641(69) at h = 0.000875 (dot-dot-dashed curve). 
Ay = 0.703(77) at h = 0.0004375 (dashed curve) and Ay = 0.687(65) at h = 0.00021875 
(solid curve). In the latter case, the quadratic fit within k £ [/c20)fc43] (dot-dashed curve) 
also is very good and looks plausible. It yields Ay = 0.694(63). In fact, all these fits at 
h = 0.0004375 and h = 0.00021875 give well consistent results within the statistical error 
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0.2 0.4 0.6 0.8 k 

Figure 2: The effective longitudinal exponent Aefr depending on the wave vector magnitude 
A; at L = 512 and h = 0.000875 (empty squares), L = 512 and h = 0.0004375 (empty 
circles), L = 512 and h = 0.00021875 (solid circles), L = 384 and h = 0.00021875 (solid 
diamonds), as well as at L = 256 and h = 0.00021875 (x). The cubic fits of the L = 512 
data within k G [fc2o,/c7o] at /i = 0.000875, h = 0.0004375 and h = 0.00021875 are shown 
by dot-dot-dashed line, dashed line and solid line, respectively. The dot-dashed line 
shows the quadratic fit of the latter data set within k G [A;2o, ^43]- 
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bars. Therefore, our final combined estimate is Ay = 0.69 it 0.10, where 0.69 is the rounded 
average value of the two fits at the smallest field h = 0.00021875, and the error bars are 
roughly estimated as ±0.10 to include the statistical error, as well as the systematical 
error due to small finite-Zi effects. This systematical error is assumed to be smaller than 
the discrepancy between the estimates at h = 0.000875 and h = 0.00021875. Our estimate 
A|| = 0.69ib0.10 is clearly inconsistent with the expected standard-theoretical value Ay = 1. 
In fact, the curves in Fig. [2] deviate away from 1, and it looks very unlikely that effective 
exponent could converge to this value as an analytic function of k. 

One has to note that corrections to scaling, proposed by the GFD theory |17j . are 
represented by an expansion in powers of A;^"^-"-, /c'^^"^!! and fe^H . The existence of a small 
correction-to-scaling exponent 2 — Aj_ is very important. It makes the extrapolation of the 
'^efF(^) plots unreliable, so that the true asymptotic value could be remarkably different 
from A|| = 0.69 it 0.10. Nevertheless, the actual estimation is well justified as a test of 
consistency of the standard theory. 

We have analyzed also the effective transverse exponent. This, however, does not give 
a better numerical evidence than those already considered in |26t [27] . We note only that 
the effective exponent (in the limit lim lim ), evaluated approximately from fits of the 

/i— >0 L^oo 

Gj_(k) data within [k,4:k], apparently, has a maximum around k ^ kj with the maximum 
value Amax ~ 1.976. Thus, A^ < Amax < 2 most likely holds, in agreement with (j32]) . 

6 Conclusions 

1. Different theoretical predictions for the power-law singularities and critical expo- 
nents of the n-vector model (or model) have been considered in Sees. [21 [3] and O 
In particular, it has been discussed how to relate the possible values of the critical 
exponents, proposed by the GFD (grouping of Feynman diagrams) theory in |16] . 
to specific n-vector models. This approach has been considered as an alternative 
method to the perturbative RG treatment. 

2. A non-perturbative proof concerning corrections to scaling in the two-point correla- 
tion function of the (^^ model has been provided in Sec. [3l showing that corrections 
to scaling proposed by the GFD theory rather than those of the perturbative RG 
theory can be correct. 

3. The known fits of the experimental specific heat data very close to A-transition point 
in liquid helium have been critically reconsidered and tested in Sec. [H It turns out 
that the overall fits with the RG correction-to-scaling exponent 9 0.529 fail to 
give satisfactory results in certain test of validity. We have demonstrated also that 
these experimental data can be very well interpreted with the critical exponents of 
the GFD theory, according to ([2]), @ and dlD. 

4. Goldstone mode singularities in the 0{n) models have been discussed in Sec. [5l 
showing that recent Monte Carlo estimates are in agreement with the theoretical 
predictions of the GFD theory. A new Monte Carlo evidence has been provided, 
according to which the statement of the (old) standard theory that G'||(k) ~ k~^\^ 
with A|| = 1 holds at /i = +0 and A; — )• in three dimensions does not look plausible, if 
corrections to scaling are such as proposed by this theory. In this case our estimation 
yields An = 0.69 ±0.10. 
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